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Abstract 

We have constructed an approximate microscopic model for the neutron dynamic structure fac- 
tor of solid methane in phase II. We expect our model to apply for neutron energies below lel^ at 
pressures near 1 bar and temperatures below 20K, where solid methane exists in a phase, called 
phase II, in which the tetrahedral hydrogens exist in a mixed phase of free rotors and hindered ro- 
tors. Certain matrix elements needed for the dynamic structure factor S(Q^uj) are adapted from a 
calculation by Ozaki et a/.[lj[2j. The model treats the effects of molecular translations, intramolec- 
ular vibrations and the free and hindered rotational degrees of freedom as uncoupled ones. Total 
scattering cross sections calculated from the model agree with the cross section measurements of 
Grieger[3] and Whittemore[4j for the incident neutron energies of 0.5 meV leV. The effective 
density of states extracted from the model was compared with the Harker & Brugger frequency 
spectrum extracted from neutron scattering measurements conducted at 22K[5j. 
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I. INTRODUCTION 



The theory of phase space compression of neutron ensembles, usually called moderation 
in neutron physics, has been well-developed for nuclear engineering and related applications 
in the energy regime which neutrons possess upon liberation from nuclei (MeV) to the 
kinetic energies of atoms in matter at room temperature ( 25 meV). Since these energies are 
large compared to the kinetic energy and the binding energy of the atoms in the moderating 
medium, the total cross section is given to a good approximation by the incoherent sum of 
the scattering cross sections from the individual atoms in the medium. In this limit, the 
theory of neutron moderation needs only to apply energy and momentum conservation to a 
sequence of collisions of neutrons of mass m and initial energy Ei with target nuclei of mass 
M and zero kinetic energy. Such an analysis forms the core of neutron moderation theory 
in nuclear reactors. The large number of collisions typically needed to reduce the energy of 
a neutron by several orders of magnitude leads inevitably to neutron emission sources which 
are incoherent. 

For neutron energies at or below the energies of motion of the moderating medium, 
however, the situation is qualitatively different. In this case, the neutron wavelength is 
larger than the separation between atoms in the medium, and in this case the scattering 
exhibits strong interference effects which depend on the details of the structure and modes 
of motion of the medium. Furthermore, the inelastic fraction of the total cross section which 
is, of course, the essential component for phase space compression, decreases as the neutron 
energy drops below the binding energy and the meV excitation energies of collective modes 
in the medium. Interaction of the neutrons with external fields are of little help in phase 
space compression: the gravitational interaction is conservative, the weak interaction is too 
weak, and the small magnetic moment of the neutron leads to inelastic effects from neutron 
spin flips in external magnetic fields in the fieV range, which is too small to be practical as 
a moderation technique for neutrons with kinetic energies in the meV range. Therefore, in 
practice, the only efficient way at present to increase neutron phase space density is through 
collisions in a medium. 

The amount of phase space compression per collision for slow neutrons in a moderating 
medium, therefore, tends to decrease as the energy decreases. Although one can in prin- 
ciple increase the phase space density in a neutron moderator indefinitely in the absence 
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of absorption by cooling the medium to T = and allowing the neutrons to undergo an 
arbitrary number of collisions, in practice the finite neutron absorption cross section of all 
media other than ^He places an upper bound on the number of collisions that can be toler- 
ated, and then once again the amount of phase space compression per collision depends on 
the microscopic properties of the medium even at T = 0. In addition the inelastic modes 
available in a condensed medium tend to freeze as T ^ 0, thereby reducing the efficiency of 
the moderating medium even further. 

The rate of phase space compression per collision is proportional to the neutron dynamic 
structure factor S{Q^uj) in the medium, where uo — Ei — Ef^ Q — ki — kf and E"^, /c/, 
ki are the final and initial neutron energies and momenta. In the theory of neutron phase 
space compression in the cold neutron regime, we are especially interested in the behavior of 
S{Q^ cj) as cj ^ 0. This is because the width of the distribution of neutrons in phase space 
that we are interested in cooling is typically large compared to the range in (Q^ou) space 
over which S{Q^ u) is large. When the latter range is defined by the dispersion relation of 
a well-defined elementary excitation of the medium such as a phonon, magnon, libron, etc., 
iS'(Q, cl;) is large only over a narrow ridge in (Q, cj) space. If the neutrons of interest for cooling 
also possessed a phase space distribution of comparable width, the obvious optimal solution 
to neutron cooling would be to match the elementary excitation spectrum to the neutron 
phase space to be cooled and reduce the neutron energy to zero by the creation of a single 
elementary excitation, with the refrigerator that maintains the moderating medium at T = 
providing the dissipation required for phase space compression. This strategy is used in so- 
called superthermal neutron moderators optimized for the creation of ultracold neutrons. 
For cooling a broader phase space distribution of neutrons to energies in the lOOfieV ^ meV 
range, however, multiple collisions are essential, and S{Q^uu) must, therefore, possess some 
strength near — > so that most neutrons in the distribution have energy losses oui in some 
sequence of collisions S{Q^ uoi) that can allow them to approach Ef 0. 

These considerations guide the search for improved neutron moderating media into some 
obvious directions. First of all, it is clear that normal phonon excitation becomes 
inefficient at low energies, since the density of states and therefore S{Q^uj) vanishes as 
u ^ 0. The abnormally large scattering cross section of hydrogen, typically more an order 
of magnitude larger than for other nuclei, makes hydrogenous materials the obvious choice 
if the number of collisions for phase space cooling is not limited by the 0.3 barn neutron 
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absorption cross section of hydrogen at 25meV and its 1/VE increase at lower energies. 

Solid methane is a special case of particular importance for cold neutron moderation 
theory for a number of reasons. First of all, solid CH4 possesses a high number density of 
hydrogen atoms. Second, the high symmetry of the CH4 molecule is ultimately responsible 
for one of its most interesting properties: namely, that a significant fraction of the tetra- 
hedron of H atoms in CH4 are free to rotate even close to T = in the solid phase. This 
property is shared by few other molecular solids such as solid hydrogen and 7-picolene [6j. 
The resulting inelastic free rotor modes are responsible for the observation, confirmed em- 
pirically through measurements of neutron spectral intensity at spallation neutron sources, 
that solid methane is the brightest cold neutron moderating medium[7j. Third, the motions 
of the rotational, translational, and vibrational modes of CH4 are sufficiently decoupled and 
the molecule is simple enough in structure that one can foresee the possibility of perform- 
ing a theoretical calculation of the dynamic structure factor S{Q^uj) of reasonable accuracy 
from first principles [3j. Such an analytical treatment in the theory of neutron scattering 
can be more easily used to gain physical insight into the cold neutron moderation process. 
In particular, one can analyze quantitatively some suggestions that have been made in the 
past to modify the CH4 medium to increase its cold neutron brightness. For example, in 
the T < 20 phase of solid methane, phase II, only 1/4 of the tetrahedral are free to rotate, 
and 3/4 of the sites undergo libration and tunneling motions. Many have speculated that 
one could improve the cold neutron brightness from a methane moderator if one could mod- 
ify the CH4 in some way so as to free up the rest of the modes. However, the librational 
and tunneling motions also possess inelastic components, and it is not obvious which of the 
two modes presenting in phase II solid CH4 is, in fact, more effective in neutron cooling, 
although the commonly-accepted folklore is that the rotational mode is the essential com- 
ponent. Fortunately the librational and tunneling energy spectra for phase II solid CH4 
have been measured [8] [9], and so the theory need not attempt to perform the delicate cal- 
culation to predict these effects from first principles but rather can simply take them from 
experiment. Yet another theoretical simplification comes from the identity of the 4 hydro- 
gen atoms, which correlates the form of the rotational and spin wave functions through the 
overall symmetry they must possess under identical particle interchange. 

Another consideration which makes this work timely is the development of new types of 
bright cold neutron sources. At high-power neutron spallation neutron sources, the radiation 
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damage in solid methane from fast neutrons and gammas is so severe that a sequence of 
events is initiated which leads in short order to the catastrophic warm-up of the moderator. 
The Low Energy Neutron Source (LENS) at Indiana University Cyclotron Facility (lUCF) 
is a long pulsed accelerator based neutron source based on (p, n) reaction in the Be target. 
The low heat load on the LENS moderator due to the greatly reduced number and energy of 
fast neutrons and capture gammas in this type of source make it possible to operate our solid 
methane moderator, and also other delicate moderating media which would be destroyed 
if operated in a MW-power spallation source environment, at temperatures well below 20K 
with infrequent warm-ups to release stored energy. In addition to its use as a cold neutron 
moderator at LENS, solid methane is planned to be used as a premoderator for an ultracold 
neutron (UCN) source at the NC State PULSAR reactor. 

In this paper, we present an approximate model for S{Q^uj) for phase II solid CH4. The 
rest of the paper is organized as follows. Section II discusses the structure of solid methane 
in phase II. Section III discusses a neutron scattering model of solid methane. The results 
including neutron cross section and spectral intensity are on section IV. We conclude and 
suggest further work in section V. 

II. SOLID METHANE IN PHASE II 

Solid methane at low temperatures has two distinct solid phases. The detailed phase 
diagram is shown in Fig. [T] At equilibrium vapor pressure, phase I appears from the triple 
point temperature (90. 6K) to the transition temperature (20.4 K) into phase II. In phase 
I, methane has a face-centered cubic (fee) structure with four orientationally disordered 
molecules per unit cell. This rotational motion of tetrahedral methane molecules remains 
slightly retarded even near the melting point of the crvstal[10j. 

The phase change is correlated with a structural shift from the point symmetric group Oh 
in phase I to the symmetrical group and Oh in phase II. In phase II, below 20. 4K, the 
center of mass of methane molecules is located on the site of the face-centered cubic lattice. 
However, a structural transition to a partial orient ational order occurs. This ordering of 
molecular orientation leads methane to possess two different sublattices, and Oh in 
phase II [lOj. 

Fig. [2] shows the partially ordered structure of phase II. It is called an eight-sublattice 
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FIG. 1: The phase diagram of sohd methane [Uj 

antiferorotational structure. The unit cell of phase II contains 32 molecules. Three-fourths 
of the molecules are hindered rotors which librate in the minima of a strong orientational 
potential of symmetry, the other 1/4 of the molecules are free rotors. 

This structure was predicted by James and Keenan on the basis of electrostatic octopole- 
octopole inter action [12]. The structure of phase II was experimentally determined by W. 
Press [13j through neutron diffraction on CD 4. In the solid CD^^ phase II is stable in the 
temperature range 22K<T<27K. Although the direct structure of CH4 was not possible 
to determine because of the high incoherent cross section of protons, they assumed that 
the experimentally measured symmetry of CD4 implies the same orientational structure of 

A. Rotational Excitation of Free Rotor: Oh symmetry 

In phase II, CH4 is in a partially ordered state with respect to the molecular orientations. 
The eight sublattice structure consists of six ordered sites with point symmetry 42m and 
two disordered sites with point symmetry 432. The disorder is caused by the partial or 
complete loss of the three dimensional periodicity of physical properties. In a molecular 
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FIG. 2: The structure of solid methane in phase II[3J. The circle one represents almost free rotation 
of the molecule. The tetrahedra denotes the orient at ionally ordered molecules in the structure. 

solid, orientational disorder can be produced if a part of the anisotropically interacting 
molecules is replaced by a isotropic unit^l4j. Meanwhile, the ordering is caused by the 
octopole-octopole interaction. The site of 432 symmetry feels only a weak crystalline field 
and the methane molecules are therefore almost free quantum rotors. The free rotor is 
aflPected more by disorder because the relaxation of orientations around a defect removes the 
compensation of octopole-octopole interaction. 

The Hamiltonian of the Extended James-Keenan (EJK) model for molecule i in a crystal 
field Vc^ from the Eq. 2.15 in Yamamoto.et a/.[llj, is 

l>4 

The numerical values of /?4 and /?6 are listed in Table. I of Yamamoto.ei a/.[llj. Ki is the 
operator of the Rotational kinetic energy. The parameter fc = 1.25 was determined by 
fitting the strength of the crystalline field [llj. The rotational wave function is given by 
the standard Wigner functions Dl^^^uj). The rotational states of a tetrahedral rotor with 
potential parameters determined for CH4 are shown in column B of Fig. [3} The energy 
levels are classified through the symmetry f F of the corresponding wave function. The label 
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FIG. 3: The energy levels [T5] of two different spin states of CH4 molecules: (a) Oh is orient at ionally 
disordered molecules (b) D2d is orient at ionally ordered molecules 

rr denotes for the 25 irreducible representation of the direct product group O x O which 
describes the symmetry with reference to the crystal-fixed frame (F) and with reference to 
the molecular fixed frame (F). The total wave function of a CH4 molecule is the product of 
a rotational wave function and a spin wave function due to the Pauli permutation principle, 
in which the total wave function must be anti-symmetric with a permutation of two protons. 
Therefore, each state is related to a representation A, T and E corresponding to the total 
nuclear spin 1 = 2 (ortho), / = l(meta) and / = O(para), respectively. The state EE cannot 
be reached by neutron interaction from the ground state AA for the reason of nuclear spin 
conservation [6j. 

The more simple description with integer concerning the energy level which begins with 
J = for the ground state has been notated by Ozaki et at.[l\. As may be seen in column A 
and B of Fig. [sj the J = 1 level is lowered from the free rotor value by 25% from crystalline 
fields. There is some ambiguity in J = 2 state concerning the energy level splitting into 
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additional states. The energy transitions between rotational energy levels are 1.09 meV 
inJ = O^J=l transition, 1.56 meV inJ=l^J = 2 transition and 2.65 meV in 
J = ^ J = 2 transition. 

B. The rotational tunneling in molecular field: D2d symmetry 

The molecule on the orientationally ordered lattice point of CH4 librates around the min- 
ima of the orientational potential. There is also the additional possibility to tunnel between 
each of the indistinguishable equilibria in two or three dimensional rotational motions of 
CH4 molecule. The rotational tunneling describes the phenomenon of the librational state 
of a molecule whose rotating atoms are indistinguishable ^6j. The splitting between the 
states is called Rotational Tunneling^^ . 

The calculation of eigenstates for using free rotational wave functions was performed 
in the framework of EJK model by Yamamoto et a/.[TT]. These procedures required a 
diagonalization of an infinite dimensional Hamiltonian matrix. Although the procedure 
converged quickly because of the weak orientational potential, the wave functions of these 
tunneling levels are only slightly different from those of the completely free rotor. 

An alternative approach was developed by Hiiller and KroU using a "pocket state" to treat 
rotational and tunneling excitation in molecular crystal[8j [9j. A systematic representation 
is shown in column C of Fig. [3) A rotational potential with n minima on two or three 
dimensional unit sphere was considered . These minima have a symmetric relation to each 
other. A set of states \i > centered at each minimum is taken as the basis set. |i > is 
the state of molecule called "pocket state" which can find itself in this basis set. Applying 
operator from the symmetry group T of the CH4 molecule, one obtains the states \n >= 
Rn\i >. Here Rn denotes one of the 12 symmetry elements of the tetrahedral group T (T=23: 
subgroup of proper rotation of point group 43m). Hiiller explicitly gave the transformation 
which diagonalizes the 12 Hamiltonian matrix < n^\H\n > of the system. For tetrahedral 
symmetry, there is just one 120° overlapping matrix element /i = /ii = /12 = /^a = /i4, as 
well as 180° overlapping matrix element H = Hi = H2 = = H4. Diagonalization yields 
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states with A^T and E symmetry with the following eigenvalues 



Ea 



D + 3H + 8h 



(singlet) 



(2) 



Et = D-H 



(three triplets) 



Ee 



D + 3H-4h 



(doublet), 



where D =< n^\H\n > is the diagonal element of the Hamiltonian matrix. In a strong 
orient ational field (in case H « h)^ the tunneling of double axes can be neglected 0. 

The oscillation frequency, called the tunneling frequency, depends on the overlap of pocket 
states which is determined primarily by the strength of the rotational state. From the Pauli 
principle, each rotational state of symmetry A, T and E is associated with a definite total 
spin of the four protons, / = 2,1,0 respectively. The energy levels of each state are 0, 
0.16 and 0.24 meV^ respectively. The librational state starts at Euij = 6.5 meV. Transitions 
between levels of diflFerent symmetry cannot be induced by phonon interaction alone, but are 
mediated by the weak dipole-dipole interaction among protons. As a consequence, the spin 
system shows slow thermal relaxation after a sudden change of the lattice temperature[17j. 

III. NEUTRON CROSS SECTION MODEL OF SOLID METHANE 

A methane molecule in a crystal shows three diflFerent degrees of freedom in the solid 
phase. The rigid H4 tetrahedron rotates about the central carbon atom. The whole molecule 
oscillates in the crystal lattice. Finally, the protons vibrate around their equilibrium posi- 
tions. 

The total neutron scattering cross section was evaluated from the scattering function 



where b denotes the scattering amplitude and ki and kf the momentum of neutron before 
and after the scattering. 

The scattering function S{Q^ cu) contains the complete dynamics of the methane molecule; 
rotation about the center, transition of the center of mass and molecular vibration. If there 
is no coupling among these modes, S{Q^ u) may be written as a convolution of the dynamic 
structure factors for the different degrees of freedom of the motion: 




(3) 



S{Q,Lj) = Srot{Q,<^) (S) Strans{Q,<^) ® Syib{Q,Uj). 



(4) 
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By integrating S{Q^Lju) over the solid angle Q and the final neutron energy the total 
cross section atot is obtained. 

If each degree of freedom is active in different neutron energies, S{Q^uj) in Eq. [3] can 
be further approximated by three different terms depending on the incident neutron energy 

En- 

For En < 6.5meV^ the protons in CH4 behave like bound particles. Only the low- 
est few rotational energy levels are thermally populated, while the molecules are in their 
translational and vibrational ground states [3j. Therefore, for the translational and vibra- 
tional motions, S{Q^Lju) can be approximated with Debye- Waller forms ey^p{—^transQ'^) and 
exp(— 7^^5(5^) respectively. The 7 stands for a model parameter 7 = j^ot + Itrans + 7vib — 
2.15 X 10-Sm2[3j. 

Then the scattering function S{Q^uj) in this energy region can be approximated as 

S{Q,UJ) = S'^ot((5,6j)exp {-{jtrans+lvib)Q^) • (5) 

For 6.5 meV < < 165 meV ^ the tetrahedron reaches its rotational "free" state 
and the librational degree of freedom is active. The cross section for rotations decreases, 
and the cross section for molecule oscillations in the crystal lattice appears , 

S{Q,UJ) = exp {--fvibQ^) Srot{Q,(^) ® Strans{Q,(^)- (6) 

When neutrons reach the energy En = 165meV^ the protons start to behave like individ- 
ual free particles. The neutron scattering is then dominated not only by multiple inelastic 
scattering, but also by vibration of the tetrahedrons. The rotational degree of freedom is 
now treated as a simple Debye- Waller form , 

S{Q,UJ) = exp {--frotQ^) Strans{Q,^) ® S^ii,{Q,Uj). (7) 



A. Rotational Excitation 

The rotational excitation of solid methane in phase II was taken from Ozaki's work rT][2] 
which has a more complete description of rotational excitation of methane. The double 
differential cross section for rotational motion can be divided into two parts: 

(Pa _ / dPa \ / (fa 

dVtdE ~ [dVtdEJ . . \dVtdE 

\ / spin— ma \ 
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Each term of Eq. [8] represents coherent and incoherent scattering in the rotational system 
as 



/ 72 \ ^ 

= hlHfp^{T)9^9f + l)jf (kr) X F^Sihu + Ef - E,) 

V / spin—ind ^ /_q 

[-jr^] = bL^P^iT)g^9f J^i^l + l)jf (kr) X G^Sihu + Ef- E,), 

V / svin—dev ^ j—n 



spin—dep ^ i—q 

in which the coherent scattering corresponds to spin-independent part and incoherent part 
to spin-dependent one. Pi{T) are the occupation probabihties for the initial states in thermal 
equilibrium 

PiyJ- ) 



'^expi-Ei/kBT) 

i 

From Eq. |9} the scattering functions of molecular rotation in each rotational state are 

oo 

SZliQ, ^kf = P^iT)9^9f + (kr) x F^^ Sihu + Ef - E,) 

5;-(g, Lo\j = MT)g,9f + (kr) x G^Sihu; + Ej - E,). 



1=0 



Here gi and gf denote the degeneracy of the initial and final states and F/'-^ and G]'^ are 
transition integral coefficient between the initial and final states which were calculated in 
Ozaki's work [l]p]. i/(kr) is the spherical Bessel function. The r stands for C — H distance 
within a methane molecule (r = 1.093 A). 

After summing over all initial and final states, the complete expression for the scattering 
function for molecular rotations in solid methane including coherent and incoherent parts is 

Srot{Q,co) = J2(^Zt{Q,^kf + sZiQ^^kf)- (11) 

B. Intermolecular Vibration : Multi-phonon Excitation 

The incoherent scattering from a Bravais lattice due to one-phonon processes is 

''^ 6LET^^('^-±S/)|Q^Vrx^<!^±|±ile-- (12) 



q,s 



where Q denotes momentum transfer and V is the polarization vector for each phonon 
mode [18]. 
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If |Q • in Eqjl2 is replaced by |Q • VT where the bar denotes that an average over a 



constant frequency surface has been taken, the incoherent cross-section may be written as 



where z{fi) is a frequency of zth phonon and cubic symmetry has been assumed. 

Since phonons of different wave vectors are independent in the incoherent approximation, 
the generahzation of one phonon scattering to the multi-phonon case is straightforward. We 
considered a number of different phonon modes changing their population by ±1 and thus 



n-phonon cross section contains a production of n-f actors of type |Q • V| |(/)| . It can thus 
be shown that the cross-section for neutron downscattering is 

SM£ = ''-^^sn(L^^— ^IFrj'^' + (14) 

Using the Debye approximation, it is possible to carry out the necessary integrations to 
obtain the cross-section cr^^^. The total inelastic incoherent cross-section is then given by 



oo 

mc 



^inelastic 



n=l 

The scattering function Strans{Q->^) is now written as 

oo 



(16) 

^-1 .--1 V-'-OO Jl r^k.T ^ I 



n=l i=\ \^-^ e^i' 



C. Intramolecular Vibrations 

From the simple harmonic oscillator model [19], the scattering function for intramolecular 
vibration is 

1 

S{Q, uj) = — - / dt exp{-iujt) x (exp(-iQ • f )exp {iQ • r(t)}) . (17) 

To calculate the response function for our model, we need to calculate the correlation 
function, 

(exp(-zQ • r)exp {iQ • f (t)}) = ]J (exp(-zQ^ • rc)exp {iQ^ • fa(i)}) 

= exp{-((Q.r)) + ((Q.r){Q.r(t)})}. 
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For an isotropic harmonic oscillator, 



((Q.r){Q.r(t)}) = ^^^^"^^^"°(;;+^^f^ (19) 



On setting t = 0, we find on the right side of the Eq. [T9| 

^ ■ = (2^) (^^^°^) = 

Using the identity, 

00 

exp(?/)cosh(x) = exp{nx)In{y) (21) 

— oo 

where n is an integer and In{y) = In{—y) is a Bessel function of the first kind. Using 
X = uj^iit + \ fil3) and v = ^v/n— a , the correlation function becomes 

exp((Q-r){Q-r(t)}) =^/„(l/)exp|na;o(«t+^^/9)|. (22) 

— OO ^ ^ 

and we obtain the scattering function as 

^ oo 

S{Q, uj) = exp(-2Vr(g) + -hul3) In{y)S{nu - nhuoo) (23) 

— OO 

where uuo is normal vibration frequency. A methane molecule has four normal vibrations 
with characteristic energies of huo = 168, 172, 376 and 391meV [ED]. 

IV. RESULT AND DISCUSSION 

A. Scattering Function S{Q,uj) of Solid Methane 

The total scattering function is the convolution of scattering functions of the three differ- 
ent modes. As mentioned, we treat this total scattering function as the linear combination 
of each scattering function in Eq. [3j Eq. |4]and Eq. [5] depending on incident neutron energy. 
This is an approximation to calculate the total scattering function when each different scat- 
tering mode dominates in different neutron energy regimes. We show the 3D plot of S{Q^ u) 
in 20K and 4K as a function of momentum transfer (Q = ^ 10 A ^) and energy transfer 
{huj = ~ 20meV) in Fig. 4(a) and Fig. |4(b)| In these figures, the two peaks in lower 



energy transfer (0 ^ hcu ^ lOmeV) represent the densities of energy and momentum trans- 
fers from the rotational modes. The long knoll around hu ~ 12meV denotes the density 
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(a)20K 




due to the multi-phonon excitation . For lower temperature, the overaU density for energy 
transfer from the rotational mode is increased. The energy transfer to phonons is essentially 
unchanged. 



B. Generalized Frequency Spectrum 

It has been remarked that the important part of the scattering law for neutron thermal- 
ization problem is the self-term defined as Van-Hove self-correlation function 



Gs{r,r) = ^ J^dr6{r + rn{0) -r)6{r - r^(r)) 



(24) 



and the scattering function have 

S{Q.^) = ^J e'^'^"^^Gs{r,T)drdT. (25) 



Eq. [24] is the velocity correlation function for the motion of an atom in the medium. The 
physical definition of this function is as follows: one images that the velocity of an atom 
is measured at a certain time ti. Then the atom is followed as it moves about the system 
and its velocity is remeasured at a later time ^2- Finally, the product of these two velocity 
measurements is taken and averaged over all values of the velocity which may be generated 
by thermal excitation of the system [2 Ij, 



'u;'expl^-^)s{Q,u)\ _ exp (- 



2kBT 



j 27r 

-P{u), 



(26) 



2M 

which defines the function P{oj)- The velocity correlation function is a complex function 
and its real and imaginary parts are related to the function P{u!) as 

Im {vq{0)vq{t)) = — J ijp{u!)sm{u)T)du 

^ / / ^ / ^^ W f°° P(cu)cOs(cUt) , 

2M Jq co8h{riuj/2kBT) 

where p{(jj) = P{(jj) ^^^^[^/^l^^r^^ is the generalized frequency spectrum and it has a condition 
to be p{uj)duj = 1. The function p{uj) can be applied in many cases and it is the same as 
the frequency distribution in a solid. But in the case of any other material, it depends on the 
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types of motion which are going on in the system. In sohd methane, the frequency spectrum 
represents aU characteristics of the rotational, translational and vibrational behaviors of a 
molecule [22j. 
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FIG. 5: The frequency spectrum of solid methane in the 20K and 4K temperatures. The Harker 
& Brugger spectrum was in the 22K 

Fig. [5] shows the frequency spectra of solid methane in two phase II temperatures, 20K 
and 4K, compared with the 22K Harker & Brugger's frequency spectrum[5j. The black 
solid line is the Harker & Brugger's frequency spectrum in 22K which was generated by 
interpolation from data. The red dotted line is the frequency spectrum at 20K and the blue 
dashed line is at 4K. The first peak around 2 meV in the frequency spectrum in 4K and 20K 
represents the linear superposition of the energy transitions from ground and the first excited 
rotational energy states to higher states in free rotational mode and the tunneling among 
ground state multiplets (three lowest states) in hindered rotational mode. The second peak 
denotes the linear superposition of the energy transition to the higher rotational states in 
free rotational mode and the librational states in the hindered rotational mode. The final 
peak around 12.9 meV is a multi-phonon mode corresponding to the Debye energy {ksTD)- 
The multi-phonon scattering and higher excitations are suppressed at low temperature, but 
the transitions probabilities from the free rotor state and ground state multiplets in hindered 
rotor increase. 



22K Harker & Brugger 

20KShin 

4KShin 
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C. Total Neutron Scattering Cross Section 



The total neutron scattering cross sections have not only been calculated from the double- 
differential cross section in Eq. [3| but also generated from the frequency spectra employing 
LEAPR module of NJOY code. 

In the low neutron energy limit, the scattering is mainly elastic scattering from free and 
hindered rotors and from tunneling among multiplet of rotational ground states in hindered 
rotor. We include the elastic coherent and incoherent scattering of hydrogen from Ozaki 
[I][2] and Grieger[3]. 

Around = lOmeV^ the total cross-section approaches 2006 in two different temper- 
atures. For En < lOmeV^ atot mainly depends on translational excitation. In the higher 
energy region, the vibrational excitation becomes more important than the rotational and 
translational scattering. 

Fig. |6(a)| and |6(b) show the total neutron scattering cross section in 20K and 4K. In 
20K, the theoretical calculation and the cross section generated from the frequency spectra 
through LEAPR module are compared to 20K measurement data [4J . The frequency spec- 



trum was used directly in multi-phonon expansion in the LEAPR module. In Fig. 6(a) 



the black solid line is the cross section calculated from the model and the red dotted line 
and blue dashed line are the cross-sections at 4K and 20K which NJOY generated from the 
frequency spectrum. Although there is a certain difference between these two cross-sections, 
the overall cross-section of model and NJOY agree well with experimental data in 20K. In 
the 4K cross-section plot, however, there is a difference between theoretical model calcula- 
tion and NJOY one. This is primarily due to the difference of elastic scattering cross-section. 
From the measurements and the theoretical calculation, it was clear that as the tempera- 
ture is lowered, the elastic cross-section a^i should be increased and inelastic cross-section 
ainei is decreased. However, the elastic scattering cross-section generated from NJOY wasn't 
changed even in the lower temperature which causes about 1006 of discrepancy between the 
theoretical calculation and NJOY across O.lmeV neutron energy in the 4K. 
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FIG. 6: The total neutron scattering cross section of solid methane in 20K and 4K. The NJOY curve 
is from the cross-section generated from the LEAPR module in NJOY through the multi-phonon 
expansion method 

D. Decoupling of the Rotational Modes 



Despite of many speculations, as we mentioned in Chapter 1, there has been no clear 
explanation on which one, between free and hindered rotations, is more effective in improving 
the cold neutron neutron flux. To approach to the answer, we assume that the one of the 
rotational modes could be turned off so that solid methane molecule has either free or 
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FIG. 7: The Elastic Scattering Cross Section in 4K and 20K from the LEAPR module and 4K 
from the model calculation. 

hindered mode only in its rotational degree of freedom in our model. In addition, the multi- 
phonon scattering and vibrational mode are still in present. From the assumption, the 
scattering function, frequency spectrum and neutron scattering cross section in each case 
are compared . 



Fig. [81 [Q] and [T0| show the scattering function, frequency spectrum and neutron scattering 
cross section of each cases in 4K system temperature. In the S{Q^ uj) plots, the higher density 
of energy transfer are mainly concentrated in low energy area for free only rotational mode. 
The hindered rotation mode shows two separate distributions corresponding the rotational 
tunneling and the librations in high energy region. These features are consistent in the 
frequency spectrum plot. The free rotation mode has about twice higher density of low 
energy transfer than hindered mode in Q ^ limit. The density of librations in hindered 
rotations are also low. This higher density of inelastic scattering of free rotation made shows 



better neutron inelastic scattering cross section in Fig(TOj Even though there has been no 
report to get all free up rotational modes of solid methane in phase II, it is clear from our 
studies that the free rotational mode of solid methane in phase II certainly show higher 
probability than hindered rotation mode in improving the cold neutron flux. 
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(a) Free Rotationl and Phonon Modes 




(b)Hindered Rotation and Phonon Modes 
FIG. 8: The 3D plots of scattering function S{Q,u) in 4K. 
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FIG. 9: The frequency spectrum of solid methane in the 4K temperatures. The Harker & Brugger 
spectrum was in the 22K 
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FIG. 10: The Scattering Cross Section of sohd methane in the 4K temperatures. 



V. CONCLUSIONS 

Theoretical studies of the neutron scattering cross section model of solid methane were 
performed. Based on the model, scattering function S{Q^Lju) and frequency spectrum p{uj) 
for solid methane in 20K and 4K are evaluated. The LEAPR module of NJOY code was also 
used for the evaluation of total scattering cross section from the frequency spectrum in the 
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20K and 4K. The total scattering cross sections are shown to be consistently explained in our 
model. The role of each rotational mode was estimated in the scattering function, frequency 
spectrum and scattering cross section by switching off the other rotational mode. A neutron 
energy spectrum measurements and the MCNP simulations will be made to validate our 
scattering model. 
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